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The study of particle motion on spherical surfaces is relevant to adsorption 
on buckyballs and other solid particles. This paper reports results for the 
binding energy of such dimers, consisting of two light particles (He atoms or 
hydrogen molecules) constrained to move on a spherical surface. The binding 
energy reaches a particularly large value when the radius of the sphere is 
about 3/4 of the particles ' diameter. 



1. INTRODUCTION 

The problem of reduced dimensionality is one of the main themes in 
condensed matter physics. One reason for an explosion of interest in this field 
is the discovery of carbon nanotubes, which are quasi-one dimensional (ID). 
A similar growth of research activity has focused on clusters, which are zero- 
dimensional from the standpoint of statistical mechanics (because no length 
diverges). A typical problem in these fields is the relation between properties 
in the reduced dimension and those in our 3D world. For example, the critical 
temperature of the liquid-vapor transition Tc increases monotonically from 
zero to a value of order the pair interaction energy as D increases from 1 to 3 
or more. One might wonder, furthermore, about the evolution of behavior for 
systems of variable effective dimensionality De// , such as the fluid adsorbed 
on a cylinder of varying radius. For very large (small) radius, Dg// = 2 (1). 

In this paper, we evaluate the binding energy of two particles, whose 
nuclei are constrained to move on a spherical surface of radius R. This is 
related to recent work exploring the evolution of films adsorbed on a Cgo 
buckyball and similar substrates It is also related to an analogous study 
of dimer binding energy on a cylindrical surface!^. That study found that 
dimers are very strongly bound if the cylinder's diameter is of order the 
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hard-core diameter a of the interacting particles . In that situation, the 
binding energy greatly exceeds that found in the 2D limit of infinite radius. 
The reason for this large binding is the existence of a favorable phase space 
condition at a particular ratio of R/cj. We report here similar behavior in 
the spherical case; that is, a particularly large binding energy occurs if 
an analogous geometrical condition is satisfied. 

While this paper focuses on for the dimer, experience with the anal- 
ogous cylindrical surface problem suggests that properties of the many-body 
system will also exhibit extrema at a particular intermediate value of the 
ratio R/(7.EE1 



2. GROUND STATE ENERGIES AND WAVE FUNCTIONS 

One may think of this geometry by imagining that one particle resides 
at a fixed position, the pole, ri =(x,y,z)=(0,0,R) while the other particle 
moves across the spherical surface in the presence of their interaction. The 
Schrodinger equation uses the Hamiltonian: 



H 



2^ 



Vl + V{r). (1) 



Here /i is the interacting particles' reduced mass, V2 is the 2D Laplacian 
(variable polar and azimuthal angles, 6 and (j)) and V(r) is taken as the 
usual Lennard-Jones(LJ) interaction with well-depth e. The coordinate r 
is the 3D separation: r=2R sin(^/2). In all cases studied here, we have 
assumed the LJ parameters: for the He isotopes, o"=2.556 A and e =10.22 
K, while for H2-H2 the parameters are a = 3.04 A and e = 34.3 K. Then, 
Eq. (1) can be reduced to a dimensionless form; the resulting equation 
provides an implicit relation between the reduced binding energy E^ =EB/e 
and the ratio R/cr. This reduced equation involves a dimensionless quantity 
h/ a^fie)^^"^ , related to the de Boer quantum parameter, discussed below; 
when this quantity is large, the reduced zero-point energy is large and the 
binding energy is small. 

A useful function is the specific area a(r), defined as the area of the 
sphere lying within a distance interval [r, r-|-dr] from ri, divided by dr. A 
simple calculation yields the following result: a(r) = 27rr0(2i? — r), where 
0(x) is the Heaviside unit step function. This result means that the most 
surface area lies at distance from ri nearly equal to 2R, but there exists no 
surface at separation r>2R. To take optimal advantage of the attraction, 
therefore, the potential minimum should lie at distance somewhat less than 
2R. One can prove the following consequence of this geometry: the classical 
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Fig. 1. Binding energies for the various He dimer problems as a function 
of sphere radius. The curves labeled "var" are variational results while the 
others are exact. 

second virial coefficient of a fluid on this surface Bciass(T) has a maximum 
when R=(t/2 and decreases monotonically with R if R> a/ 2. One might 
investigate whether similar behavior occurs in the quantum case. Then, the 
low temperature behavior is dominated by the spectrum of bound states, if 
any are present. 

We have determined the ground state energy (-Es) for this problem in 
two ways, variational and exact, for a number of different systems, including 
^He -^He, ^He -^He and ^He -^He, as well as H2-H2 . An interesting question 
is this: what ratio of R/cj gives the largest binding energy E^. Since no 
attraction is provided by the LJ potential unless R/a >l/2, one might guess 
that the largest binding occurs when R/a Furthermore, we expect 

the system to attain the largest reduced energy Eg when the particles are 
most classical; a completely classical system would have Eg =1 since the 
particles would maintain a spacing equal to the minimum in the potential. 
The criterion for "most classical" is the largest value of the de Boer quantum 

parameter: A* = -j^— , where m is the mass of the interacting particles, 

assumed to be the same. The values are A*= 3.1, 2.7 and 1.7 for the sequence 
^He, ^He and H2, respectively . 

Figure 1 presents the exact and variational energies for the various He 
cases. There is seen to be a maximum binding energy of order 1 to 2 K, with 
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Fig. 2 

Fig. 2. Normalized probability density as a function of polar angle for the 
various cases studied here. 

specific values 1.97 K, 1.38 K and 0.88 K for the sequence ^He - '^He, ^He - 
^He and ^He - '^He, respectively. This trend manifests the expected decrease 
of binding energy with increasing A*; a similar behavior was reported by 
Kilic et al. for dimers confined to either flat surfaces or to the interior of an 
spherical cavity The maximum binding occurs with R/cj 0.72, 0.76 
and 0.80. This is analogous to the case of dimer binding of He isotopes on 
a cylindrical surfacd^, with a (somewhat smaller) maximum binding at R/a 
=0.65. The smaller value of R/cj in the cylindrical case is expected because 
the a(r) cutoff of long range contributions to the attraction, mentioned above, 
is not present there. 

Figure 1 also shows the result of a variational calculation, using a trial 
wave function (r) of the form: 



The two variational parameters b and n were optimized for each radius. 
The resulting energy in the ^He - ^He case, shown in Fig. 1, has b=1.80 and 
n=2.50 at the optimal radius R/o"=0.80. Note that the variational result for 
the binding energy is significantly lower than the exact binding energy, while 
the optimized radius is nearly identical to the exact one, R/o" =0.80. An 
analogous situation was found in the problem of the dimer on a cylindeJ^. 
There, the optimized variational binding energy was 0.084 K, while the exact 
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energy was 0.218 K. Interestingly, the two calculations yielded the same 
optimal radius in that case too. 

Figure 2 compares the probability densities for the various cases studied 
here. The characteristic feature is an increasing localization near the south 
pole as the system becomes more classical, as implied by the a(r) discus- 
sion above. Note that the probabilities derived from the variational wave 
functions are more diffuse than the exact solutions. 

In Fig. 3 we present the exact binding energy for the H2 - H2 dimer. In 
this case, the most classical of the group reported here, the maximum Eg is 
0.5. This value is much closer to the classical value (1) than that for ^He - 
^He (E^=0.2). One might not have anticipated such a large increase in 
since the A* parameters are not so different ( 2.7 vs 1.8). Note, however, 
that for the He isotopes, just a 15% difference between values of A* leads to 
a factor of two difference in E^. Evidently, the quantum parameter plays 
a critical role in this problem because the hard core potential confines the 
atoms to a fairly small domain on the surface. 



3. SUMMARY 



Wc have found particularly large binding energies for dimers of He and 
H2 when the particles are constrained to move on spherical surfaces having 
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R/cr ~3/4. This behavior is quahtatively similar to that found previously 
for cylindrical surfaces. A possible next step would be to explore the quan- 
tum virial expansion on this surface. Such a theory would take a different 
form from the phase shift approach used for extended systems. Prom the 
experimental point of view, the present work suggests the intriguing possi- 
bility of spectroscopic observations of these dimer states, e.g. in a mixture 
of Ceo molecules and He atoms. A quantitative prediction would require a 
study of both the excited state spectrum and the role of the third dimension, 
associated with the radial motion, which is ignored in the present work. 
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